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This article completes a former work where part of the 0{p^') low-energy constants entering in 
the TTTT scattering were estimated. Some resonance contributions were missed in former calculations 
and slight differences appeared with respect to our outcome. Here, we provide the full results for all 
the contributing 0{p^) couplings. We also perform a reanalysis of the hadronic inputs used for the 
estimation (resonance masses, widths...). Their reliability was checked together with the impact of 
the input uncertainties on the determinations of the chiral couplings and the scattering lengths Oj. 
Our outcome is found in agreement with former works though with slightly larger errors. However, 
the effect in the final values of the aj is negligible after combining them with the other uncertainties. 
Based on this consistency, we conclude that the previous scattering length determinations seem to 
be rather solid and reliable, with the 0{p^) low-energy constants quite under control. Nevertheless, 
the uncertainties found in the present work point out the limitation on further improvements unless 
the precision of the couplings is properly increased. 

I. INTRODUCTION 

In a previous work [|], ||, we provided a set of predictions for some of the 0{p^) low energy constants (LECs) 
Ti related to the TTTr-scattering amplitude The 0{p^) chiral perturbation theory couplings r2, ... rg were 

determined there in the large- A^c" limit Q by means of once-subtracted partial- wave dispersion relations, where they 
were provided in terms of the ratios of widths over masses, Vn/M^ and Vj^/M^. The limit of a large number of 
colours Nc — > oo ||^, is a key ingredient of the study, becoming the strong dynamics greatly simplified and being 
the dominant contribution provided by the tree- level meson exchanges. At large Nc, the relevant resonances for 
TTTT-scattering are the 1=1 vector and the uu -f dd component of the scalar with 1=0. We will denote these large- A^c 
states, respectively, as p and a all along the paper. We will also consider just the contribution from the lightest 
multiplets: the single resonance approximation will be assumed. 

The new predictions in Ref. Q found that some scalar resonance contributions to the 0{p^) LECs had been actually 
missed in former estimates This produced small variations on the TTTr-scattering lengths of the order of the 
current errors However, in order to make a thorough analysis, we complete our former study and provide the 
remaining 0{p^) LECs contributing to TTTr-scattering. In the standard 0{p^) chiral perturbation theory calculation 
the scattering lengths Oj depend on ri,... rg. Alternatively, if they are determined through the dispersive method in 
Ref. one needs to provide r^^ instead of r^ and rg. 

In addition, we have decided to redo the whole analysis of the resonance inputs. This has enabled us with at least 
a minimal control of the uncertainties in the LECs and the corresponding scattering lengths. This allows a better 
understanding of what are the relevant parameters that determine the low energy scattering and the main sources of 
errors. The vector mesons are found to fit very well within a U (3) large- A'^c multiplet and their properties are quite 
under control. On the other hand, the current knowledge of the lightest large-iVc scalars is rather poor. Thus, our 
revised analysis of the TTTr-scattering lengths is also motivated by the need of improving the current picture of the 
lowest lying hadronic resonances [0, g], ^ • 
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In Section II, we present the former LEG determinations and the reanalysis that we propose. The phenomenological 
inputs required for the determination of the LEGs are studied in full detail in Section III, keeping a careful control of 
the possible uncertainties. Based on this, the values of the 0{p^) low-energy constants are estimated in Section IV, 
leading to a series of new predictions for the tttt scattering lengths in Section V. The various results and conclusions 
are gathered in Section VI. Finally, some more technical details are relegated to the Appendices. 



II. LEG ESTIMATES THROUGH RESONANGE SATURATIONS 

The current calculations in chiral perturbation theory (xPT) have already reached the two loop order |3[ Q . How- 
ever, in order to extract quantities such as the scattering lengths one needs to eventually input the 0{p^) low-energy 
constants. Since in many cases it is not possible to determine these couplings directly from the phenomenology, one 
needs to extract their values through alternative procedures. One of the most usual ones is to consider large- A^c" esti- 
mates based on phenomcnological chiral lagrangians |l^, 0. However, the control that one has on these lagrangians 
and their large-TVc estimation may be unclear. Thus, the uncertainties of the observablcs under consideration might 
be larger than what is actually quoted. 

In this work, we propose the comparison of three sets of estimates of the low-energy constants. First we will review 
the currently used values (set A) ^ and then we will present the new numbers after taking into account the scalar 
meson contributions that were missing in former works (set B) ||^. However, for sake of consistency, we will also redo 
the analysis of the hadronical inputs and provide the newly calculated LEGs (set G). 

• Set A: 



This is the group of estimates commonly employed in nowadays calculations ||^, U]. The xPT couplings are 
assumed to be determined by the resonance exchanges provided by the phenomenological lagrangian 



^ {u^,u^' + x+) 



+ ^( V^5V^5) - iil/K SS ) + Cd{ Su^u^ ) + c,n{Sx+ ) 



(1) 



where (...) stands for trace in flavour space, S and account respectively for the scalar and vector multiplets. 
The tensor w'^ contains the chiral pseudo-Goldstone and x± i^i in addition, proportional to the light quark 
masses. Their precise definitions can be found in Refs. |, |, |l§. Notice that since the lagrangian docs not 
contain resonance mass-splitting terms, the masses M/j coincide with their chiral limit values, which we denote 
as Mr- Nevertheless, it can be straightforwardly extended to include the splitting due to the quark masses. 
This will be studied with full generality in the next sections. From the comparison of the decays p — > tttt, 
K* Ktt and other processes, Ref. [|j provided the set of parameters 



Mv = 770 MeV , gv = 0.09 , = -0.03 , 

Ms = 983 MeV , c„ = 42 MeV , Cd = 32 MeV . 

Taking these inputs and the phenomcnological lagrangian ([|), they produced the LEG estimates 



(2) 



-0.6 X 10"^ 
-1.0 X IQ-^ 



1.3 X 10-", 
1.1 X 10-4. 



-1.7 X 10"4, 
0.3 X 10"\ 



(3) 



where the authors already accounted in these numbers the contribution to the SU{2) LEGs coming from the 
kaon and eta loops H : 
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These will be also included in our analysis in order to have a clearer comparison with the outcomes obtained in 
Refs. |,|. 

In addition to this, the dispersive approach used in Ref. ||^ also required the TTTT-scalar form- factor coupling rs^, 
which was estimated to be 



-0.3 X 10" 



(5) 



with a negligible contribution r 



• Set B: 



11527r^m 



0.03 • 10 from kaon and eta loops 



K 



However, the resonance estimate from Refs. y, |^ was incomplete and some relevant contributions were actually 
missing in their calculation. The presence of the operator c,„(S'x+) in the scalar meson lagrangian produces 
a tadpole term proportional to the quark mass Although its effect vanishes in the chiral limit, its 

contribution becomes relevant at subleading chiral orders. 

The 0{p^) LECs are recalculated here in full detail, keeping all the possible contribution. No term of the 
corresponding chiral order under study is neglected. By means of the partial-wave dispersion relations proposed 
in Refs. |Q, ^, most of the 0(p^) LECs in the TTTr-scattering (r2, ... rg) are now fixed in terms of the ratios. 



Ml 



Ml, 



R 



M 



R 



M 



R 



l + c.R^+lR^,+0{ml) 



M 



l + pR^+0{m%) 



(6) 



where M n and Tn stand for the chiral limit of Mr and Tn, respectively. The constants a^,, Pr, 7i? are quark 
mass independent and rule the corrections in the ratios. The formal expressions for r2, ...rg Q are provided 
later in Sections IV A| - |rV^ . The resonance lagrangian allows us then to compute the resonance masses and 
widths at large Nc in terms of the resonance couplings. Thus, the LEG predictions can be utterly rewritten in 
terms of the latter. Using exactly the same inputs (||) of set A, one obtains then slightly different values 



= 18x IQ-"^. 



0.9 X 10" 



= -1.9 X 10" 



(7) 



remaining and rg unchanged, i.e., = rf and = as in Eq. (|3|). Although the variation in ra,... rg 
with respect to Set A is relatively mild, r2 changes by an order of magnitude. This points out the need of a 
more detailed investigation of the stability of the estimation under modifications in the inputs. 

To complete the calculation we needed to extract ri (for the direct xPT calculation [||) and rg^ (for the 
dispersive method |^). These two couplings were out of the reach of the partial- wave dispersion relations 
proposed in our former works ||l|, ||. Here we have used the lagrangian (|^). We have performed the explicit field 
theory calculation of the tttt scattering A(s, t, u) and the scalar form factor !Fs{s). The details of the calculation 
are shown in Appendix A. 

The corrected determinations derived from the phenomenological lagrangian (|l|) result 



ri = 



rS2 



9v 



Ml 



Mt 



(8) 
(9) 



which for the set A inputs in Eq. (0) yield 



r? = -2.1 X 10~^. 



-0.3 X 10" 



Our prediction agrees accidently the numerical value r^^ reported in Ref. 
found to be small. 



(10) 

and the variation in ri is also 
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• Set C: 

The analysis can be further refined. In addition to performing the full eomputation of the LECs (without 
neglecting any contribution), a more careful examination of the phenomenological inputs also seems convenient. 
A more general description of the mesonic interactions is required. For instance, the lagrangian considered in 
Eq. (|l|) does not take into account the fact that in physical QCD there is a mass splitting within the U{3) 
resonance multiplet. This quark mass effect utterly contributes to the scattering lengths at the same order as 
the terms already included in the determinations from sets A and B. A similar thing happens with the resonance 
widths, which accept a more general quark mass splitting pattern. 

Actually, it is possible to describe most of the U(3) breaking without relying in any particular resonance 
lagrangian realization. Nevertheless, without lost of generality, we will study the vector resonances in both the 
Proca and antisymmetric formalisms p^ . Although the two representations are physically equivalent, small 
discrepancies may appear when considering corrections only up to a given order. The slight difference that 
appears when rearranging the experimental information from one formalism to the other will serve us as an 
estimate of the residual error due to higher order contributions. We will also perform a general analysis of 
the splitting in the resonance masses and widths. This will allow us to fix their dominant chiral corrections, 
being refiected in a more accurate estimate of the LECs r^,... rg. However, it will be impossible for us to have 
a full control of the subdominant chiral corrections -that are going to affect ri, r2 and rs2~-- Therefore, we will 
have to utterly rely on estimates of these LECs based on phenomenological resonance lagrangians like that in 
Eq. (|l|). This Set C reanalysis of the resonance parameters will be shown in detail in the next Section. 



III. PHENOMENOLOGY OF THE RESONANCE PARAMETERS 
A. Mass splitting up to 0{m%) 

In the large- A'c' limit, the mass splitting of the resonance multiplets can be described at leading order by a single 
operator 



2 

which leads at large Nq to the mass eigenstates 



"'RR) + e^{RRx+), (11) 





^1 


= 


-Ae^ml + 


0{m%) , 




1 
2 


= Ml 


-4e^m^ + 


0{m%), 


Mj2 


2 




- 4e^. {2ml 


-ml) + 



O(m^) . (12) 
and m\%+''''^ = Mi=i. 

In the case of the vector resonance multiplet, can be fixed through the physical p{77Q) and AT* (892) masses, 
respectively Mp = 775.5 ± 0.3 MeV and Mk' = 893.6 ± 1.9 MeV §]. This provides the estimates 

Mv = 764.9 ±0.5 MeV, ejf, = -0.217 ± 0.004 . (13) 

These values and Eq. ([ij) lead to a prediction for the remaining two states, Mj'I^q = 998 MeV and jyf = ~ 

776 MeV. This is in pretty good agreement with the (/)(1020) and ^^(782) masses, respectively = 1019 MeV and 
M^j = 783 MeV. Hence, alternatively, one can think of extracting the chiral limit resonance mass and the splitting 
term by means of the p(770) and 0(1020) masses, with = 1019.455 ± 0.020 MeV @: 

My = 763.7 ±0.5 MeV, e^, = -0.2414 ± 0.0022 , (14) 

with the prediction = 906 MeV, in relatively good agreement with the experimental AT* (892) mass [| . 

The physical masses are indeed slightly different to their large- A^c- values. Thus, the relations in Eq. ([2|) are 
not exactly fulfilled for the experimental inputs as they gain subleading contributions in 1/Nc- The combination of 
Eqs. (^) and ( p^ provide our most reliable estimate of the large- A'^p parameters. Considering an interval that covers 
both determination, we will take as inputs from now on the values 

Mv = 764.3 ± 1.1 MeV, e"^ = -0.228 ± 0.015 , (15) 
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whose errors exceed those that come from purely the experimental uncertainties in Eqs. (|T^) and (p^. 

In the case of the scalars, it is relatively straightforward to identify the lightest 1=1 resonance with the ao(980): 
A'//=i = 984.7 ± 1.2 MeV §. The remaining states of the multiplct are more cumbersome to pin down. The scalar 
spectrum provided in the large- A^c study of Ref. jl^ leaves the light and broad scalar resonances a and k out of any 
classification multiplet. In order to avoid the problem of the mixing of iso-singlet scalars, the analysis is performed 
with the 1=1/2 state, whose mass in the large- A^c limit is also poorly known. The broad k(800) seems to be a 
possible candidate although the first clear 7=1/2 scalar resonance signal is provided by the i^o(1430) Hence, we 
take the conservative estimate = 1050 ± 400 MeV, which ranges from the k up to the i^o(1430) mass. This 

leads to the predictions 

Ms = 980 ±40 MeV, e^ = -0.1±0.9. (16) 



B. The splitting of tlie vector resonance decay width up to 0{m%) 

The calculation of the quark mass corrections to the width is slightly more complicate and one needs to make some 
assumption on the structure of the interaction. 

If we use the phenomenological lagrangian (|^), the vector decay width into two light pseudo-scalars, V 0i</>2, 
shows the structure 



4o \ ^2 



1 + ev— ^ 

2M,, 



H 2 



0(m%) 



(17) 



with the phase-space factor 



PV12 



1 - 



(mi + 



1 - 



(mi - m2y 



(18) 



The Fi are the physical decay constants of the pseudo-Goldstones [Ft^ ~ 92.4 MeV and Fk — 113 McV) and they 
appear due to the large- A^c wave function renormalization of the light pseudo-scalars |||, My and corresponds, 
respectively, to the physical vector and pseudo-scalar masses. Xytttt is the Vcjiiipi coupling that rules the decay 
amplitude in the chiral limit (e.g. gv in Eq. (|^) Q) whereas ey rules the quark mass corrections to the amplitude 



(provided for instance by the couphng in Eq. (y) [^). Cvi2 is the corresponding Clebsch-Gordan: CpTr-n ~ 1, 
Ck*Ktt = 3/4 and C^^-f^ = 1- The vector mass scaling n is an integer that depends on the resonance lagrangian 
under consideration: n = 5 in the Proca formalism ^ and n = 3 in the antisymmetric tensor representation |l^, ^ . 
Finally, one power of pvi2 comes from phase-space and the other two powers_^are due to the J^= 1 spin polarization 
sum of the squared modulus of the amplitude, jA^y^i; 
denote the vector meson polarizations. 

Thus, in the Proca field realization from Eq. dll) one has 



|2e^p5'7W|2 = M^plr^^\M\^, where here the 



5, 



A 



Vtttt 



gv ■ 



4V2/x/5v 



(19) 



On the other hand, in the antisymmetric tensor representation employed in Resonance Chiral Theory (R^T) |T^, 
one has Xvtttt = Gy, n = i and the quark mass corrections ey provided by a combination of some extra V4>(j) operators, 
i)^ {V^v{x+,u'^u''}) , iXg (Vfj.^u^^x+u'^ ) 0- The detailed discussion on the quark mass splitting operators in the 
antisymmetric formalism is relegated to Appendix Abusing of the notation, we will denote as the effective 
combination of resonance chiral theory couplings such that for the antisymmetric formalism [n = 3) we still keep 
ev = AV2,yXv^^ = 4.V2fJGv 

In order to remain as general as possible, both vector formalism scalings, n = 5 and n = 3, are discussed. The 
following inputs will be used all throughout this paper to estimate the uncertainties of our results fl: 



TV = 92.4 ±0.3 MeV, Fk - 

= 137.3 ±2.3 MeV, 
= 149.4 ± 1.0 MeV, Tk-- 



113.0 ± 1.0 MeV, 
= 495.6 ±2.0 MeV, 
= 50.6 ±0.9 MeV, 



(20) 



together with the mass parameters My and derived in Eq. (ha 
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The combination of the experimental K* and p widths allows us to fix the parameters \v-kit and ev'- ^ 



\ n— 5 
\ n— 3 



- gv = 0.0846 ± 0.0008 
Gv = 63.9 ± 0.6 MeV, 



,n=5 



,n=3 



0.01 ±0.09, 

0.82 ±0.10, 



for n=5 , 
for n=3 . 



For n = 5; this corresponds to the coupling = (0.2 ± 1.3) x 10 ^, to be compared to the determination 
from Ref. |0] ^. In the case of the antisymmetric formalism, one obtains f^^^ = 9.3 ±1.1 MeV. 



(21) 
(22) 

-0.03 



The only difference between the Proca (n=5) and the antisymmetric tensor formalism (n=3) is that in the 



latter the p 



TTTT width carries the factor Gy 



Instead, in the Proca case, this is replaced by 



1 



n=5 ™„ 



glMy 



1 + {el=^ 



M, 



ri=3 ™, 

. Up to higher chiral order corrections, one finds a pretty 



good agreement between both of them, getting identical results for Gv = 63.9 MeV and gvMy — 64.7 MeV, and for 



^n— 3 



0.8 and (e 



n— 5 



Ael) = 0.9. 



The obtained values of Xv-mr and ev lead us to a prediction for the (f) 

^j^j. = 4.1 ±0.8 MeV, 



KK decay width: 



(23) 



KK 



for both Proca and antisymmetric formalisms. This result is perfectly consistent with the experimental value 
3.54 ± 0.10 MeV At large-TVc, the (/>(1020) is identified with the ss component of the / = vector. 

Likewise, the determination of the coupling Xvtttt produces for the p tttt decay width the chiral limit prediction. 



Ytttt 



487rF'* 



Pp = 181.9 ±3.0 MeV, (n=5), 
Pp = 177.8 ±2.5 MeV, (n=3), 



(24) 
(25) 



where we have used the value of F = 90.8 ± 0.3 MeV derived in Sec. HID. 



C. The decay width for the scalar resonance 

In the case of scalar resonance decays, the constraints to the structure of the width are given just by an overall phase- 
space factor pvi2 and the F^^^ factors due to the wave-function renormalizations of the pseudo-Goldstones |p^ . 

The issue that arises here is the current poor knowledge on the structure of the scalar sector in the large-iVc 
limit. One needs then to be assisted by some phenomenological lagrangian. In the present work, we will assume 
that our scalar interactions are ruled by the action from Eq. (|^) together with the mass splitting operators 
from Eq. { V\\). Although many works have tried to pin down the values of the two scalar couplings, and in 
Eq. (|l|) dio], [P^, |l^, |l|, |l|, 11], stiU there are no widely accepted resuhs. 

The oo 7r°?7 decay seems to be the only reliable process to estimate Cd and c™. However, this only fixes a 
combination of the two and one needs to add extra theoretical information. Thus, the study of the high-energy 
behaviour of the Kir, Kr], Krj' scalar form-factors performed in Ref JlGf provided the constraints Cm = Cd and 
4^CdCrn = ■ However, the first constraint seems a priori less reliable, as the analysis did not include all the possible 
quark mass operators contributing at that order (these extra terms in the lagrangian were derived later in Ref. |pT[). 
On the other hand, the relation AcdCm = is more solid, as it stems from the high-energy analysis of the scalar 
form- factor in the chiral limit. 

Hence, the constraint AcdCm ~ F^ and the oq — > ttt; decay width are used to fix the values of and c^. Thus, the 
lagrangian (^ yields 



P _ / Ficos^i - \/2F8sin6>8 y eg A/^^ Pa^^,, | -, •■-tt ' , I 



FiF8Cos(0i-08) ; 24^f2 \ Ml^ ' 



(26) 



^ Indeed, the system of quadratic equations has two possible solutions for A and ev- In one case is much larger (ev — —10) than in 
the other. Since this corresponds to huge chiral corrections, we will only keep the small solution ~ 1). 

^ If, alternatively, wo consider the same inputs as Ref. pi together with Fk = llOMeV, then one gets fx/dv — —0.03, leading to 
fx ^ -0.003. 
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where the lagrangian from Eq. (|]) yields 65 + 2 = 2cm/cd- The 77-77' mixing is given in the basis of the octet r^g and 
singlet r;i by il|,i|,Pl: 



7/ \ ^ J_/ Fg cos 6*8 -Fi sin 6*1 \ f 778 
?/ y F V ^8 sin 6*8 Fi cos 61 M 771 



(27) 



We consider the inputs = (1.1 ± 0.1) F^, 6»i = (-5 ± 1)°, Fg = (1.3 ± 0.1) F^, 6*8 = (-20 ± 2)° [gl|, |2|, ||, 
Tao—Trr; = 75 ± 25 McV, TO,, = 547.9 MeV together with the inputs considered in previous sections. Most of the 
error in the next determinations comes from our poor knowledge on the ao(980) decay width. 

Relying on the estimate of the ao rjii'^ decay width from the phenomenological scalar lagrangian (^ |l^ and 
the theoretical constraint ^CdCm = F"^ |l^, Q , one gets the values 

Cd = 26±7MeV, = 80 ± 21 MeV . (28) 

As the error is dominated by the large ao(980) width uncertainty, we find that the determination of Cd and Cm is 
completely unaffected by whether one considers the constraint AcdCm = F^ or the approximation AcdCm = F^. The 
chiral limit F of the pion decay constant F-^ is computed in the next section. Indeed, the determination of Cd is rather 
model independent since the {es + 2) term in Eq. ( |26| ) has little impact in the decay width. Thus, if the (es + 2) 
term is neglected one gets = 31 McV. 

Taking this into account, the a tttt decay width is given at large Nc by 



3 Pa-nn cl A 
IGttF^ 




(29) 



where by a we denote the iso-singlet (uu + dd) scalar, without strange quark content. The quark mass correction eg 
and the chiral limit of the width are then provided by 



Cm. , \ . , „ ^ 3 AI c 



es = 2 — -1=4±3, = — 2 — f = 600 ±300 MeV. (30) 

Cd J IdnF'^ 



D. Chiral corrections to 

In order to extract the LECs related with to^ corrections, we will need to know the quark mass dependence of F-^, 
which can be parametrized in the general form 



F^ 



F 



1 + SF, 



(2) 



SF, 



(4) 



Oiml) 



(31) 



The pion decay constant F^ appears in the calculation when one takes into account the pion wave-function renormal- 
ization Z^r that occurs at large Nc, which obeys F-^ = F ■ Z^^^"^ [|[ Q 



The scalar lagrangian in Eq. (Ill) ||10| and the mass splitting from Eq. (|12|) yield 



5F, 



(2) 



F2 



5F, 



%CdCm ( 3CdC, 



(4) 



F2 



F' 



F2 



\<OCdCmCr, 
F2 



(32) 



Although (5F(2) has the most general structure, this is not true for (5F(4). If a more general set of scalar operators 
A'^Os, X^^Oss, X^^^Osss were allowed in the lagrangian ||ll[, (5-F(4) would gain a whole series of new contributions. 
Substituting our former determincitions of c^^ Cm — F^/icd and e^^ in Eq. (|3|), one gets 



SF, 



(2) 



SF, 



(4) 



-5±5, 



(33) 



where the large uncertainty comes both from Cm and the 25% error in Cd- For (5F(4) one can use indistinctly Cm = 
F'^ /Acd or Cm = F^ I^Cd, as the difference results negligible and goes to the next order in the to^ expansion. By means 
of Eqs. (|3l] ) and (|3^) it is then possible to recover the large-iVc value for the pion decay constant in the chiral limit, 
F = 90.8 ± 0.3 MeV, in agreement with former large- A^c xPT determination pi|. 
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E. NLO chiral symmetry breaking parameters or, Pa 

Combining the information obtained from the mass and width sphttings, one can now extract the corresponding 
quark mass corrections in the ratios Vr/M^ and Tr/M^ defined in Eq. (||): 

M ^ 

av = 2ev - 2{n - 3) - A5F(^) =t " 6 , 

Ms 

M ^ 

I3v = 2ev - 2(n - 5) - ASF^^^ =T " 6 , 



Ms 



as = 265 - 4(5F(2) 
Substituting the same experimental inputs as before one obtains 



(is = 265 + 4ef„ - 4<5F(2) - 2. (34) 



ay = -7.5 ±0.3, /3y = -8.4±0.3, for 71 = 5, 

ay = -6.8 ±0.3, /3y = -7.7±0.3, for n = 3. (35) 

and for the scalar 

as = 2±7, /35 = 2±8. (36) 
Although an is not needed in the present 0{p^) LEG study, it is provided here for sake of completeness. 

F. NNLO chiral correction 'yn 

The parameter -appearing in Tj^/M'^ at 0(m^)- is even more complicate to determine than an and P^. 

At next-to-next-to- leading order (NNLO), the resonance masses may suffer corrections due to the chiral operators 

-^{RRx+x+) + -^{RX+RX+), (37) 
2M^ 2Mj, 

which combined with the leading and NLO operators in Eq. ( [Til ) yield the pattern 



= a4 - 4e,^."4 - 4(e«,i+?«,2)=^ , (38) 



MU/2 ^ - Ae'^vii, - =2 C,i (2"iK - 2mf^"4 + O + e^,2 < , 

Mr 

Mf:l ' ^mI- 4e«(2m^-m^) - 4(ef„,, + e^,^) ^^^4^^ , 

Mr 

and j\,f 1^^+'''') = Af/=i. Actually, the NNLO multiplet splitting will be only relevant for the analysis of the vector 
parameter in the Proca formalism (n = 5). In the remaining cases, = i + 2 ^ill ^'^^ appear. Since 
these NNLO quark mass corrections may enter in serious competition with those NLO in 1/Nc, an analysis of the 
vector spectrum in order to fix seems unreliable. Thus, we consider just some conservative bounds. It is not really 
possible to obtain a tight constraint from the / = 1 or / = i. The most stringent bound comes from the / = (ss) 
state: 

— 2 

K\ < o 2 Ie,nl ^ 0.3, (39) 

2m^ - 



where we have demanded that the NNLO could not overcome the NLO contribution to the (/)(1020) mass (in the quark 
mass expansion). 



9 



Likewise, the expansion of the p and a resonance widths up to NNLO in is given by the parameters tR: 

n 2 



- p — ^TTTT 



Vtvtt 



1 + + 0{ml) 



M 



V 



M 



V 



M 



(40) 



A way out to estimate these chiral corrections is the phenomenological lagrangian of Eq. (n]) 



icm{cd - c,„) My 



16cl^{cd -Cm) , 8(c„ 



f 4e^ 



CdF' 



Cd 



(41) 



with ey = 4:\/2f^/ Xvtt-k and the mass sphtting operators ( pi] ) also taken into account. The fuU unexpanded expression 
for the widths can be found in Ref.O]. Using the inputs of former sections, one obtains 



?y = -0.0 ±0.3 (for n=5), 
es = -30 ±40. 



ev = -2.8 ± 1.7 (for n=3). 



(42) 



Gathering all the different contributions together, it is now possible to get the NNLO correction to the ratio 



iv 



2?y - 2(n- 3)e^ 



(4): 



M 



+ey - 4(n - 3)e;,,ey - 8eySF^2)=^ - 12ev^ + 2{n^ - 8n + 15) 



-10(<5F(2))- 



4 



— 2 — 2 

12(n - 5)el + 2ASF(^2) =t + 8(71 - 3)e^,5F(2) =^ 



(43) 



7s 



2^5 - 4<5F( 



(4) 



+e| - 4es - 8esSF(^2) + 8SF(^2) + W{6F^2)f - Sei - 2 . 



(44) 



Notice that the O(m^) corrections to ((5f(4)), the resonance masses (e^) and widths (ej?) are presented in the 
first line of Eqs. (^)-(Q). The results of Eqs. (|43|)-(|4^) arc actually model dependent since there exist several other 
resonance operators that may contribute to the NNLO parameters (5i^(4), ey, es ||ll|. Nevertheless, if one relies in the 
model of Eq. (nl) the inputs considered in former sections lead to the chiral symmetry breaking terms 



7v = 30± 10 (for n=5), 
75 = -50 ±70. 



7y = 19±8 (forn=3). 



(45) 



A more detailed observation of the NNLO parameters (first line of Eqs. (^)~(^)) shows that the impact of the NNLO 
mass splitting term is negligible. On the other hand, (5-F'(4) and es are responsible of most of the error in 7s. In 
the vector case, the effect of the width splitting term ey is not dominant, whereas the Ft^ NNLO correction (5f(4) is 
responsible of roughly the 75% of the total uncertainty in jy. The large errors we have in our input numbers for (^i^(4) 
and es reassure us in the validity of our result, even if the uncertainty of the jr covers a rather conservative interval. 
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IV. LOW-ENERGY CONSTANT DETERMINATION AT 0{p^) 
A. The determination of rs, rg 

The couplings , rg only depend on the chiral limit values of the resonance masses and decay widths |^ and, 
therefore, are the most reliably determined low-energy constants: 

r5 = — =7 \ _7 , (46) 

3Ms My 

re = _7 • (47) 

My 

B. Extraction of rs, r4 

In the case of the couplings and r4, all one needs are the resonance widths and masses in the chiral limit and 
the first correction to the ratio Fr/A/p, this is, /^r 

15. X - / IVly 



_ ■ _ -. (i + r^)^ (48) 



The coupling depends only on the vector resonance. Likewise, the scalar contribution is quite suppressed in due 
to the large numerical coefficient in front of the vector term. 

C. Extraction of r2 

Compared to the determination of r^, ■■■r(,, the coupling r2 carries larger theoretical uncertainties since the NNLO 
parameters 7s and 7^ enter into play. The prediction for r2 derived from partial-wave dispersion relations is given 
by [| 

r, = 2r, + f 1 + ^ + + ^ (7584+ 1248;3v + 1447v) , 

Mg \ 3 6 J My 

(50) 

where rp provides F^^ at NNLO in Q. It is related to our previously defined (5-F'(4) through rp ~ 2£3^4 + 
F'^5F(^4) / M g . The phcnomcnological lagrangian (|l|) produces the value of <^-F'(4) in Eq. (p2|), which in combination 
with 2^3^ = 32c2„crf(c,„ - Cd)/Ms yields 

rp — 4 ^ =4 ■ lO-Lj 

Mg Mg 

Substituting our former phcnomcnological inputs, one gets rp = {—1 zt 3) ■ 10~^. 

D. Extraction of ri and rs2 

Within the framework of tttt partial-wave sum-rules proposed in Refs. ||l|, ^, it is not possible to make any prediction 
for T"!, as they arc based on once-subtracted dispersion relation and this coupling produces just a constant contribution 
to the TTTT-scattering amplitude. 
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The couplings ri (from TTTr-scattering) and rs2 (from the tttt scalar form-factor) must be extracted through alter- 
native procedures. The values from Eqs. (||) and (H) were based on the lagrangian (^ but without accounting for the 
mass splitting effect. If the latter is included, our phenomenological model produces for ri the new predictions 



12c^ 



32(cd . 



m3 



l + ev + 



Mc 



16cdC,„(8c^-17cdC„H-12c2j 32(cd - c™)^^^ 



16G|^ 



1 + ey 



2e: 



(52) 



The detailed calculation is relegated to Appendix A 1 



Likewise, after taking into account the mass splitting in the scalar form- factor calculation (Appendix [A 2| ), 
becomes 



rs2 



8c^{c^-Cd)F^ _ 32cgc^ 



16cdCm-F^ s 

4 

Mo 



(53) 



E. Saturation scale uncertainty 

One last problem to face is the fact that the large- A^c estimate of the LECs does not carry any renormalization 
scale dependence. However, it is possible to find a so called "saturation scale" fis such that the physical LEG r[(/i) 
agrees numerically with r^'^^°° for = Hs- 

Since the standard comparison scale is /io = 770 MeV, the coupling r\{^) must be run from /^s up to /io- The 
possible difference between these two scales introduces an uncertainty in our determination. The way considered here 
to account for this lack of knowledge is to observe the variation for a wide range of /i, which in other works is usually 
taken to be in the range 500-1000 MeV |, |, |l^, H: 

Ari''^ = 3 X 10""* , Ara''-' = 4 x 10""* , Arg''-' = 3 x 10"^ , 
Ar4''= = 0.05 x 10"'' , Arg^^ = 0.5 x 10^"* , Are^' ^ 0.05 x 10"'' , 

Ars2''= = 1.5 X 10"'', (54) 

where the running is completely fixed in xPT by the expressions given in Appendix B in terms of F and the 0{p'^) 
invariants li = -0.4 ± 0.6, £2 = 4.3 ± 0.1, h = 4.4 ± 0.2 @ and £3 = 2.9 ± 2.4 |26|. 

Sometimes, it is argued that there must exist a common saturation scale for all the LECs, both 0{p^) and 0{p^). 
This scale /is is very often identified with the rho mass, /is — Mp. However, explicit one-loop calculations with 
resonance lagrangians show that, though the loops typically produce logarithms of the form In or ln||r-, the 
combinations of them appearing for each LEG are not necessarily the same [|7[ 29, Hence, although in 



general terms /is ^ Mb,, nothing ensures that the saturation scales must be exactly identical, nor that their value 
must be equal to Mp Thus, the uncertainty of every LEG is computed and accounted separately in the 

present work. 



F. Summary of low-energy constants 

The values for the low-energy constants f[(/i) and the corresponding errors are gathered in Table |. We choose 
the standard comparison scale /i = 770 MeV and include the uncertainty from the saturation scale estimated in the 
former section as an error. 

Some remarks about the present predictions arc in order. The LECs and rg are the most model-independent 
ones since they are determined by just the chiral limit of the resonances masses and widths (see Eqs.(B^) and (p7[)). 
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ND est. 


set A 


set C 


set C 




PI 


m 


(n=5) 


(n=3) 


10* • rl 


±80 


-0.6 


-14± 17±3 


-20± 17±3 


10* • rl 


±40 


1.3 


22 ±16 ±4 


7± 10±4 


10* • rl 


±20 


-1.7 


-3± 1±3 


-4 ± 1 ± 3 


10* • rl 


±3 


-1.0 


-0.22 ±0.13 ±0.05 


0.13 ±0.13 ±0.05 


10* • rl 


±6 


1.1 


0.9 ±0.1 ±0.5 


0.9 ±0.1 ±0.5 


10* • rl 


±2 


0.3 


0.25 ±0.01 ±0.05 


0.25 ±0.01 ±0.05 


10* -rS, 


±1 


-0.3 


1 ±4± 1 


1 ±4± 1 



TABLE I: Different predictions for the 0{p^) LECs r[(^) for /i = 770 MeV: The first column presents the order of magnitude 
estimate based on naive dimensional analysis [^2|; In the set A column we show former estimates from Refs. ^; in the last 
two columns, one can find the values for the present reanalysis, both in the Proca (n=5) and antisymmetric vector formalism 
(n=3). The first error derives from the inputs and the second from the uncertainty in the saturation scale. 

Our determinations for r^, and rg are consistent with those in Ref. The constants r^ and r4 depend on the 0{m'^) 
corrections in the ratio Ffl/Afl , which are still rather under control. All our predictions are in agreement within 
errors with those reported in Ref. ||3|] except for the small deviation found for r4, mainly due to the slight discrepancies 
in the coupling mentioned in previous sections. All this points out the little model dependence of our LEG estimate. 

On the other hand, r2 is partly determined by the NNLO corrections 77? to the ratio Vr/M^. This makes it 
one of the least controlled LECs. Roughly half of the total error of the LEG is due to the NNLO corrections (5-F'(4) 
and £5 in 75. Thus, the uncertainty of r2 is found to be dominated by the scalar mass parameters, which is due in 
our case to our poor knowledge of the 7 = 1/2 scalar mass. 

In general, the value of all the couplings that can be determined by means of the partial wave dispersion relations 
(r2,... rg) 0, ^ is found to be dominated by the vector contributions. Since this sector is rather under control, this 
reassures us in the robustness of our calculation. 

Finally, the couplings ri and r^^ need to be estimated through a phenomenological model. In the case of the 
lagrangian the error of ri is mainly due to our ignorance on the scalar resonances couplings Cd and c,„. This is 
directly originated in the large uncertainty of the ag tttj decay width. The vector contribution is much smaller and 
essentially negligible. The error in stems mainly from our poor knowledge of the ao(980) width and the scalar 
mass splitting (which derives here from the large uncertainty of the 1 = 1/2 scalar resonance). 

In summary, we find in general a good agreement with former determinations and at the same time we are able to 
provide a reliable estimate of the error. This will help us to establish the relevance of the 0{p^) LEG contributions 
to the scattering lengths in next section. 
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Set A § 

(xlO^^) 


Set C (n=5) 
(xlO--') 


Set C (n=3) 
(xlO^^) 


&8 

lOag 
lOfcg 
10a\ 
Wbl 


1.2 

4.1 
-3.7 

14 
-0.3 

0.7 


2.3 ±1.3 ±1.7 

4± 1 ±4 
-4±4± 1 
-13±8±2 

3.4 ±2.2 ±0.6 
2.6 ±0.5 ±1.5 


-0.3 ±1.0 ±1.7 
2±1±4 
-4±4± 1 
-3±5±2 
1.7 ± 1.3 ±0.6 
3.7 ± 0.6 ± 1.5 



TABLE II: Contribution of the 0{p^) LECs r[(/i) to the scattering lengths Uj and effective ranges bj in the direct xPT 
approach |^] for /i = 770 MeV. The first column (set A) shows the results from Ref. with their same inputs Ftt = 93.2 MeV, 
m,r ~ 139.47 MeV. The last two columns show the predictions based on our phenomenological reanalysis (set C) for the 
Proca (n=5) and antisymmetric tensor formalisms (n=3). The first error derives from the inputs and the second one from the 
saturation scale uncertainty. 



V. SCATTERING LENGTHS 
A. Direct xPT calculation 

The TTTT-scattering lengths were first calculated in xPT up to 0{p^) in Ref. |^. At order m^, the contribution from 
the r[(/i) LECs was found to be 





ml 
327rF6 


[5r[ 


+ 12rS 


+ 48r^' 


+ 32r2 + 192rg 




ml 


r2 + 


Url + 


12r^± 


72r^ - Srg] , 


2 1 


ml 
167rF6 


K- 


h 16r^] 








ml 
SnF^ ^ 




±24r^ 


- IGrg] 






ml 
247rF6 


V2- 


[-8r^ + 






bWr, 


ml 
6nF^ ^ 




±3r^± 







In order to work with dimensionless quantities, we have multiplied the results in Ref. |^ for fefj, 6g, a\ by to^, and for 
b\ by m^. Our estimate of the contributions to the scattering lengths and effective ranges for the standard comparison 
scale /J, = 770 MeV is given in Table both for the Proca (n=5) and the antisymmetric vector formalism (n=3). In 
order to have a clearer comparison of this quantities, the outcome from Ref. |^ is also provided. The first error in 
Tabic |l| comes from the phenomenological inputs and the second one from the saturation scale uncertainty in the 
0{p^) LECs, given in Eq.(|5|). 

An important part of the subdominant quark mass corrections has been pinned down quite accurately through the 
comparison of the decays of the different resonances in the multiplet. However, there exist a series of new operators 
(e.g. the vector resonance operator Ag' in Ref. ||ll|) whose couplings cannot be extracted in an independent way. If 
present, they could enter in effective combinations gvi'triTvY^ and Cd{m^Y^] these would be what we would really 
determine and denote respectively as gv and Cd, being used later instead of them in the LEG computation. A more 
detailed discussion is relegated to Appendix |C| In any case, the study of these operators remains beyond the scope of 
this article. 



B. CGL dispersive method 

In Ref. 1^] , Colangelo et al. combined the NNLO chiral perturbation theory computation of the scattering lengths |^ 
with a phenomenological dispersive representation. This allowed them to produce one of the most precise dctermina- 
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tions of the scattering lengths. These were expressed in terms of some dispersive integrals, the pion quadratic scalar 
radius the 0{p^) coupling £3 and a set of 0{p^) LECs (ri, ... r4, rs^)- These last contributions are actually the 



most poorly known and the aim of this paper. Following the work of Ref.| 
lengths that depends on the inputs r[(/i): 



I, we extract the part of their scattering 



7m ^ 



TO, 



327ri^^ 
ml 



[5r 



28rl - Urs,] , 



4r.c. + 4r4 



(56) 



where the Co 



can be extracted from the Cj provided in the Appendix C of Ref. Q . The LEG appears in this 
analysis |^ because the scalar radius (?'^ )5 is incorporated as an experimental information in order to fix the 0{p^) 
constant £4. On the other hand, the couplings and rg disappear here with respect to the standard xPT analysis 
of Ref. j3|, being their information encoded and replaced by the different dispersive integrals. It is easy to realize 
that the resonance contributions ajl^ carry the same ri and r2 dependence in both the direct calculation in 

Eq. ( ^5|) Q and the dispersive study in Eq. ( ^6|) |6|. The comparison of this reanalysis and the predictions from 
Ref. [pL is shown in Table [II for the standard reference scale /i = 770 MeV. 





Total: Ref. |] 




Set C (n=5) 


Set C (n=3) 




(xlO-^) 


(xlO~^) 


(xlO^^) 


(xlO-^) 




220 ±5 


0.0 ± 1.0 


1.0 ± 1.5 ± 1.0 


-1.6 ± 1.5 ± 1.0 




-444 ± 10 


0.4 ±2.0 


0±4±2 


0±4±2 



TABLE III: The first and second columns show, respectively, the total scattering lengths and the ri contribution to them in the 
dispersive method from Colangelo et al. where the authors used the ri in Eq. (2.3), Ftt ~ 92.4 MeV and — 139.57 MeV. 
The last two columns show the reanalyzed quantities (set C) for the Proca (n=5) and antisymmetric formalisms (n=3) for the 
usual scale /i — 770 MeV. There, the first error derives from the inputs and the second one from the saturation scale uncertainty. 



We find that the largest contributions to the Oq and Oq errors are produced in similar terms by ri, r2, ^3 and rg^ 
On the other hand, the impact of r^ on both the value and uncertainty of the scattering lengths results negligible. 



VI. SUMMARY AND CONCLUSION 

This article concludes the previous work from Ref. The analysis of the uncertainties of the 0{p^) LECs 
contributing to the tttt scattering lengths ^ is completed here. Former estimates based on the large- A"c limit and 
resonance saturation have been revised. Nevertheless, a series of uncertainties escape to the control of the analysis 
carried in the present article. Most of the computations are rather model independent as they are based on general 
resonance properties; however, at some points we had to rely on phenomenological lagrangians |l5j. 

All this allowed the estimate of the 0{p^) LECs. For the standard comparison scale ^ = 770 MeV, one obtains 
after combining the results from Proca and antisymmetric formalism the values 

rj; = (-17±20) x 10"*, = (17±21) X 10-'', = (-4 ± 4) x 10"* , 

= (0.0±0.3) X 10-^ = (0.9±0.5) X 10-^ rg = (0.25 ± 0.05) x 10"'' , 

r^s, = (1±4) X 10"^ (57) 

The combination of the results from the Proca (n=5) and antisymmetric formalism (n=3) yield for the dispersive 
method ^ the final resonance contribution, 

lO^aglr, = 0±3, ^O^allr^ = 0±5. (58) 

This is in perfect agreement with the former estimate although the detailed analysis of the phenomenological 
inputs casts a slightly larger uncertainty. 

Following the analysis of global uncertainties of the scattering lengths in Colangelo et a/.'s 1^, we have verified that 
the global uncertainties for Oq and Uq are not largely modified. Thus, the replacement of the values from Eq. ( ^8|) in 
the total scattering lengths leads to the updated predictions 



a^o = 0.220 ±0.005, 



-0.444 ±0.011, 



(59) 
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where the total uncertainties and central values are essentially unchanged with respect to the previous determinations 
= 0.220 ± 0.005 and lOag = -0.444 ± 0.010 [|. 

This work provides a reliable and solid estimation of this part of the 0{p^) calculation. It sets clear limits to the 
size of the 0{p^) LEG contributions, slightly conservative in some occasions. A special attention has been put on the 
quantification of errors and their precise source. It reassures us in the reliability and precision of the current scattering 
length determination based on the dispersive approach of Ref. Nonetheless, although the global uncertainties of 
the scattering lengths are barely affected by our new predictions of the 0{p^) low energy constants, the total errors 
have almost reached those from ajln- This points out the difficulty of further improvements in the accuracy of the 
TTTT scattering lengths unless the 0{p^) LEG uncertainties are conveniently reduced. 
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APPENDIX A: FULL CALCULATION BASED ON THE RESONANCE LAGRANGIAN 



1. Extraction of ri: nn scattering 

The resonance lagrangian from Eq. (|^) |l^ has been used to calculate the TTTr-scattering amplitude at large Nc- 
The vector resonances were described there by means of Proca four- vector fields V^,. In addition, we also took into 
account the vector and scalar mass splittings. The amplitude 7r+(pi)7r^(p2) ^ ^"(P3)7''"(?'4) results then 



A{s, t, u 



I Proca 



+ 



u — s 



gvt + 4V2/x 



[cdS - 2cdml + 2c„m^] ' 



(Al) 



Ft{u-M') 



gvu 



4V2/x 



with s = [pi +P2)^, t = (pi — ps)^, u [pi —piY. We denote by p and a respectively the 1=1 vector and the / = 
scalar with uu + dd content, i^^ ~ 92.4 MeV is the pion decay constants, which deviates from its chiral limit F in 
the way prescribed in Eq. ( |3l] ) due to the scalar resonance tadpole in the Cm operator of the resonance chiral theory 
lagrangian (|^) |l^ . The same happens with the physical pion mass and its value at leading order in the chiral 
expansion, m^, which are related through 



m 



2Bom 



8cm{cd - Cm) ml 



(A2) 



with m^^^^ the u/d quark mass in the isospin limit. 

Alternatively, if one employs the antisymmetric tensor formalism to describe the vector resonances ||l^, |ll|, ^ the 
scattering takes the form 



A(s, t, u) 



Antii 



Ft 



[cdS - 2cdml + 2c,„m^] ' 



F^{t-AI^^ 

en=3 



Gv + 4%/2/^=3^Va^' 



F^{u-M') 



where the meaning of Gy and /" will be discussed later in detail in Appendix O 



(A3) 
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The resonance expressions can be now compared to the xPT ampUtudc , 



niis^ rriiU — uY s'^ s{t — uY 



(A4) 



We have preferred to express everything in terms of F rather than i^^r hi order to make the chiral matching more 
transparent ||] . 

The LECs arc extracted by matching A{s, t, u)^^'^ and the chiral expansion of the resonance amphtudc A{s, t, u)^'^^ 
in Eqs. (Al) and (A_3). Since we are interested in ri, s = t = and u = 4mJ turn out to be the most convenient 
choice of momenta. The matching of the term yields for the Proca formalism (n=5) 



ri + 2rF + 16r4 - 8ii = ^ 64 



3 — =2 =T 



=74 — ~ CdCm + Cj) , (A5) 

Ms 

with the term — 128/-^i^^/il/y in the right-hand side absent in the antisymmetric tensor formalism (n=3). However, 
for the case of our phenomenological lagrangian, its numerical value is absolutely negligible. Thus, the vector con- 
tribution comes mainly from the LEG r4 = ^Y/tI^ (l + ey — (n — 5)eJ^ — i which is extracted from the 

large-TVc partial wave analysis in a very reliable way B. Likewise, although r^^ is obtained by means of the resonance 
lagrangian (P, the large errors we found for i5i^(4) in Eq. (|3^) make its estimate slightly conservative. For £3 we use, 

£3 = 4^^^^^^^^ , (A6) 
Ms 

where we converted the SU{3) large- A^'c' estimate from Ref. [|lO| into SU{2) LECs |3^. Possible pseudo-scalar 
resonances have been neglected. Putting all this information together, one obtains the value for ri reported in the 
text in Eq. (IsT 



2. Extraction of rs2- tttt scalar form factor 



We proceed now to the calculation of the pion scalar form-factor by means of the resonance chiral theory la- 
grangian (|l|) P, p^, . We repeat the calculation in detail, including also the contributions from the scalar tadpole 
operator Cm p3, |13|, and obtain 



2Bn 



O 



Ml 



(A7) 



Matching the resonance description in Eq.(A7) and the xPT result for the tttt scalar from-factor ||^ leads to the LECs 

rs,=AF^^^, (A8) 
Ms 



rs2 



C3t4 



8Cm(Cm-Cd)F^ 32c^Cd{cd - 2Cm) l6CdCmF'^ s 



M, 



M, 



(A9) 



with rs3 consistent with the value from Ref. jj]. Substituting the large- A'c value of £3 from Eq. (|A6D and £. 
AcmCd/ M s leads to the prediction for rs2 shown in the text in Eq.(p 
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APPENDIX B: 0{p^) LOW-ENERGY CONSTANT RUNNING 

The variation of the r[ (^) with ^ is given by the general pattern [|[ Q 



with the scale invariants 



(167r2)2 



In^ 



2\2 



(16^^) 



K 



,(L) 



193 104- 112- 

h + -^i2 - 6^3 + 2^4 : 



27 9 
556 68^ 



9 

248 
~9~ 



^2 



7^3 - 24 . 



if. 



755 100: 



44- 



108 



9 



K 



il - -£2 

108 9 9 



29 
432 



107- 
"36" 



(L) 



K 



AL) 



S2 



= -^-^h-'^l2. 

432 36 36 
148 62- 

— + y £1 + 2^3 + 2^4 : 



In^ — - In^ — 



i^f ^ = -20 , 



K. 



(2) 



K. 



(2) 



if. 



(2) 



if, 



(2) 



407 
~' 

232 

2 

3' 



iff' = ^ 



85 
18 



5 

6 ' 



4:' = - 



166 
""9~ 



(Bl) 



(B2) 



APPENDIX C: CHIRAL CORRECTIONS TO THE DECAY WIDTH 



In the article, the splitting in the vector decay widths due to quark mass effects was parametrized at NLO by one 
single parameter ey (or equivalently /^). This was true for the Proca resonance lagrangian in Eq. However, in 
the case of a more general hadronic action, one single parameter ey does not seem enough to describe the NLO chiral 
corrections to the decay widths. In the Proca field realization, in addition to the operator x-] ), there 

could exist higher derivative operators such as ( u^u'^} ), which also contribute to the chiral corrections to 
the decay widths at NLO in m^. Nonetheless, we will focus our digression about higher terms in the lagrangian in 
the antisymmetric tensor formalism, where the possible operators related to the NLO quark mass corrections to the 
vector width have been already constructed in Ref. |11[: 



(CI) 



However, not all these three couplings are observable in the partial decay widths we analyzed in the text, which 
now become 



PV12 Q 

48 TT fI F| ^ 



off^ 



Gcff 
V 



{ml + ml) + 0{m%) 



(C2) 



with the corresponding Clebsch-Gordan Cvi2 for the channels under consideration, CVtt = 1, Ck-Ktt = 3/4 and 
^4>KK ~ 1- Instead of the 7Ti7r-independcnt coupling Gv ■, which we had before in Eq. (|l7|), one has now the effective 
combination 



Gcff 
V = 



(C3) 



The combination of resonance couplings (2Ag' -I- A]^) rules the splitting of our studied decays and determines the chiral 



symmetry breaking parameter /. 



n=3 



(2Ar + Aro)M;. 
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Hence, the p nn, K* — > Ktt and (f> KK partial decay widths only allow the determination of G^, not its 
chiral limit value Gy- A similar situation would happen with the scalar sector and the Proca field description for 
vectors, where the presence of higher operators could make us observe in our analysis some effective combinations c^^ 
and , instead of the quark mass independent couplings and gv- In any case, the difference is assumed to be 
small in the present work and its further study remains beyond the scope of this article. 
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